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ABSTRACT 

Some possible methods to determine at the LHC masses of 
SUSY particles are discussed. 

I. INTRODUCTION 

If supersymmetry (SUSY) exists at the electroweak scale, it 
should be easy at the LHC to observe deviations from the Stan- 
dard Model (SM) such as an excess of events with multiple jets 
plus missing energy Et or with like-sign dileptons £ ± ^ ± plus 
Et [0 §, §. Determining SUSY masses is more difficult be- 
cause each SUSY event contains two missing lightest SUSY 
particles \\, and there are not enough kinematic constraints to 
determine the momenta of these. This note describes two pos- 
sible approaches to determining SUSY masses, one based on 
a generic global variable and the other based on constructing 
particular decay chains. 

The ATLAS and CMS Collaborations at the LHC are consid- 
ering five points in the minimal supergravity (SUGRA) model 
listed in Table I below [Q]. Point 4 is the comparison point ex- 
tensively discussed elsewhere in these Proceedings. For this 
point a good strategy at the LHC is to use the decays \2 ~* 
X\f- + ^ to determine the mass difference M(x°) - M(\i) [fj- 
For higher masses, e.g. Points 1-3, this decay is small, but 
xl - X ih - x\bb, xf -» X?^± - X°iQQ, and X ° 2 -> U - 
Xi^ provide alternative starting points for detailed analysis. 
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Table I: SUGRA parameters for the five LHC points. 
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H. EFFECTIVE MASS ANALYSIS 

The first step after discovering a deviation from the SM is to 
estimate the mass scale. SUSY production at the LHC is dom- 
inated by gluinos and squarks, which decay into jets plus miss- 
ing energy. The mass scale can be estimated using the effective 
mass, defined as the scalar sum of the p^'s of the four hardest 
jets and the missing transverse energy Et, 

M e s = pr,l + PT.2 + PT,3 + pT,i + Et ■ 
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Figure 1: Point 1 signal and backgrounds. Open circles: signal. 
Solid circles: ti. Triangles: W — ► hi, tv. Downward triangles: 
Z — ► vv, tt. Squares: QCD jets. Histogram: all backgrounds. 
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Figure 2: Signal and SM backgrounds for Point 2. See Fig. 1 
for symbols. 



ISAJET 7.20 [g] was used to generate samples of 10K events 
for each signal point, 50K events for each of tt, Wj with W — * 
ev, fiv, tv, and Zj with Z — > vv, tt in five bins covering 50 < 
p T < 1600 GeV, and 2500K QCD events, i.e., primary g, u, d, 
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Figure 3: Signal and SM backgrounds for Point 3. See Fig. 1 Figure 5: Signal and SM backgrounds for Point 5. See Fig. 1 



for symbols. 



for symbols. 
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Figure 4: Signal and SM backgrounds for Point 4. See Fig. 1 
for symbols. 



s, c, or b jets, in five bins covering 50 < pr < 2400 GeV. The 
detector response was simulated using a toy calorimeter with 

EMCAL 10%/ VE + 1% 

HCAL 50%/ VE + 3% 

FCAL 100%/VE + 7%, |r?| > 3 . 

Jets were found using a simple fixed-cone algorithm (GETJET) 
with R = [(Ar/) 2 + (A^) 2 ] 1 / 2 = 0.7. To suppress the SM 
background, the following cuts were made: 

• E T > 100 GeV 

• > 4 jets with p T > 50 GeV and p T ,i > 100 GeV 



• Transverse sphericity St > 0.2 

• Lepton veto 

• E T > 0.2M cff 

With these cuts and the idealized detector assumed here, the 
signal is much larger than the SM backgrounds for large M e g, 
as is illustrated in Figs. 1-5. 

The peak of the M c g mass distribution, or alternatively the 
point at which the signal and background are equal, provides a 
good first estimate of the SUSY mass scale, which is defined to 
be 

-Msusy = min(Mg, M 5r ) 

(The choice of M^ R as the typical squarkmass is arbitrary.) The 
ratio of the value M c r for which S = B to Msusy was calcu- 
lated by fitting smooth curves to the signal and background and 
is given in Table II. To see whether the approximate constancy 
of this ratio might be an accident, 100 SUGRA models were 
chosen at random with 100 < mo < 500 GeV, 100 < m 1 / 2 < 
500 GeV, -500 < A < 500 GeV, 1.8 < tan/3 < 12, and 
sgn/i = ±1 and compared to the assumed signal, Point 1. The 
light Higgs was assumed to be known, and all the comparison 



Table II: The value of M e g for which S = B compared to 
^SUSY, the lighter of the gluino and squark masses. Note that 
Point 4 is strongly influenced by the Et and jet pr cuts. 

"Point M c ff(GcV) M SVSY (GeV) Ratio" 
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1420 
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Figure 6: Scatter plot of MgrjsY = rnin(Mg, Ms) vs. M e g for 
randomly chosen SUGRA models having the same light Higgs 
mass within ±3 GeV as Point 1. 
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Figure 7: Ratio M c ff/MsusY from Fig. 6 



models were required to have Mh — 100.4 ± 3 GeV. A sample 
of IK events was generated for each point, and the peak of the 
M e g distribution was found by fitting a Gaussian near the peak. 
Figure 6 shows the resulting scatter plot of MgusY vs - M e g. 
The ratio is constant within about ±10%, as can be seen from 
Fig. 7. This error is conservative, since there considerable con- 
tribution to the scatter from the limited statistics and the rather 
crude manner in which the peak was found. 

HI. SELECTION OF h -> bb 

For Point 1 the decay chain \2 ~ ¥ Xih> h — ► bb has a large 
branching ratio, as is typical if this decay is kinematically al- 
lowed. The decay h — ► bb thus provides a handle for identifying 
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Figure 8: M(bb) for pairs of & jets for the Point 1 signal (open 
histogram) and for the sum of all backgrounds (shaded his- 
togram) after cuts described in the text. The smooth curve is 
a Gaussian plus quadratic fit to the signal. The light Higgs mass 
is 100.4 GeV. 



events containing x^'s @l- Furthermore, the gluino is heavier 
than the squarks and so decays into them. The strategy for this 
analysis is to select events in which one squark decays via 



and the other via 



xll; 



giving two b jets and exactly two additional hard jets. 

ISAJET 7.22 ^ was used to generate a sample of 100K 
events for Point 1, corresponding to about 5.6 fb -1 . Back- 
ground samples of 250K each for tt, Wj, and Zj, and 5000K 
for QCD jets were also generated, equally divided among five 
Pt bins. The background samples generally represent a small 
fraction of an LHC year. The detector response was simulated 
using the toy calorimeter described above. Jets were found us- 
ing a fixed cone algorithm with R = 0.4. The following cuts 
were imposed: 

• E T > 100 GeV 

• > 4 jets with p T > 50GeVandp Tj i > 100 GeV 

• Transverse sphericity St > 0.2 

• M efi > 800 GeV 

• E T > 0.2Af cff 

Jets were tagged as b's if they contained a B hadron with 
Pt > 5 GeV and i] < 2; no other tagging inefficiency or b 
mistagging was included. Figure 8 shows the resulting bb mass 
distributions for the signal and the sum of all SM backgrounds 
with pT,b > 25 GeV together with a Gaussian plus quadratic 
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Figure 9: The smaller of the two bbj masses for signal and back- 
ground events with 73 < M(bb) < 111 GeV in Fig. 8 and with 
exactly two additional jets j with pt > 75 GeV. The endpoint 
of this distribution should be approximately the mass difference 
between the squark and the \\, about 542 GeV. 



fit to the signal. At a luminosity of 10 33 cm~^s~ 
will have a &-jet tagging efficiency of 70% for a rejection of 
100 [[!]]. Hence, the number of events should be reduced by a 
factor of about two, but the mistagging background is proba- 
bly small compared to the real background shown. The Higgs 
mass peak is shifted downward somewhat; using a larger cone, 
R = 0.7, gives a peak which is closer to the true mass but wider. 

Events were then required to have exactly one bb pair with 
73 < M(bb) < 111 GeV and exactly two additional jets with 
Pt > 75 GeV. The invariant mass of each jet with the bb pair 
was calculated. For the desired decay chain, one of these two 
must come from the decay of a single squark, so the smaller 
of them must be less than the kinematic limit for single squark 
decay, M(u R ) - M{x\) = 542 GeV. The smaller of the two 
bbj masses is plotted in Fig. 9 for the signal and for the sum of 
all backgrounds and shows the expected edge. The SM back- 
ground shows fluctuations from the limited Monte Carlo statis- 
tics but seems to be small near the edge, at least for the ideal- 
ized detector considered here. There is some background from 
the SUSY events above the edge, presumably from other decay 
modes and/or initial state radiation. 

IV. SELECTION OF W -> qq 

Point 1 also has a large combined branching ratio for one 
gluino to decay via 



9 -> QLq, qL -> xtq, X.f XiW ± , W 



and the other via 



qq, 



giving two hard jets and two softer jets from the W. The branch- 
ing ratio for qi, — > XiQ i s small for Point 1, so the contributions 
from g — > q^q and from <?l9l pair production are suppressed. 

The same signal sample was used as in Section III, and jets 
were again found using a fixed cone algorithm with R = 0.4. 
The combinatorial background for this decay chain is much 
larger than for the previous one, so harder cuts are needed: 

• E T > 100 GeV 

• > 4 jets with p T \.2 > 200 GeV, p T 3A > 50 GeV, and 

7?3,4 < 2 

• Transverse sphericity St > 0.2 

• Afcff > 800 GeV 

• E T > 0.2M eff 

The same ^-tagging algorithm was applied to tag the third and 
fourth jets as not being b jets. Of course, this is not really feasi- 
ble; instead one should measure the 6-jet distributions and sub- 
tract them. 

The mass distribution M34 of the third and fourth highest pt 
jets with these cuts is shown in Fig. 10 for the signal and the sum 
of all backgrounds. A peak is seen a bit below the W mass with 
a fitted width surprisingly smaller than that for the h in Fig. 8, 
note that the W natural width has been neglected in the simu- 
lation of the decays. The SM background is more significant 
here than for h — ► bb. Events from this peak can be combined 
with another jet as was done for h — > bb in Fig. 9, providing an- 
other determination of the squark mass. Figure 10 also provides 
a starting point for measuring W decays separately from other 
sources of leptons such as gaugino decays into sleptons. 
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Figure 10: M34 for non-& jets in events with two 200 GeV jets 
and two 50 GeV jets for the Point 1 signal (open histogram) and 
the sum of all backgrounds (shaded histogram). 



9 -> qn.q, qn -> Xxi ■. 



V. SELECTION OF x° 2 



Point 1 has relatively light sleptons, which is genetically nec- 
essary if the Xi is to provide acceptable cold dark matter [Q], 
Hence the two-body decay 



X°2 



is kinematically allowed and competes with the X2 ~~ * Xi^ de- 
cay, producing opposite-sign, like-flavor dileptons. The largest 
SM background is ti. To suppress this and other SM back- 
grounds the following cuts were made on the same signal and 
SM background samples used in the two previous sections: 

• M cS > 800 GeV 

• E T > 0.2M e fl 

• > IR = 0.4 jet with p T ,i > 100 GeV 

• £ + £~ pair with p T ,£ > 10 GeV, r/£ < 2.5 

• I isolation cut: E T < 10 GeV in R = 0.2 

• Transverse sphericity St > 0.2 

With these cuts very little SM background survives, and the Mu 
mass distribution shown in Fig. 1 1 has an edge near 
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If Mu is near its kinematic limit, then the velocity difference 
of the £ + £~ pair and the Xi is minimized. Having both leptons 
hard requires M^/M? ~ M^jM~ t Assuming this and M^o = 

2M^o implies that the endpoint in Fig. 1 1 is equal to the Xi 
mass. An improved estimate could be made by detailed fitting 
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Figure 11: Mu for the Point 1 signal (open histogram) and the 
sum of all backgrounds (shaded histogram). 
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Figure 12: M tt ^o for events with 86 < Mu < 109 GeV in 
Fig. 11, using p^a = M^o /MuPu for the Point 1 signal (open 
histogram) and the SM background (shaded histogram). 



of all the kinematic distributions. 

M™ ax - 10 GeV < M„ < M™ ax , 



Events were selected with 
and the xi momentum was 



calculated using this crude x? mass and 

p^a = (M^JMu)pu. 

The invariant mass M U j^o of the £ + £~, the highest pt jet, and 
the xi was th en calculated and is shown in Fig. 12. A peak is 
seen near the light squark masses, 660-688 GeV. More study is 
needed, but this approach looks promising. 
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